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We present the kinematic variable, mT2i which is in some 
ways similar to the more familiar 'transverse-mass', but 
which can be used in events where two or more particles 
have escaped detection. We define this variable and describe 
the event topologies to which it applies, then present some 
of its mathematical properties. We then briefiy discuss two 
case studies which show how mT2 is vital when reconstruct- 
ing the masses of supersymmetric particles in mSUGRA-like 
and AMSB-like scenarios at the Large Hadron Collider. 



1 Introduction 

Reconstructing R-parity conserving supersymmetric events will be difficult 
at the Large Hadronic Collider (LHC) because of the following factors which 
limit our knowledge of the event: 

• two massive particles have escaped undetected, 

• the masses of these particles are unknown. 
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Invisible 
final state 



Figure 1: Schematic representation of a simple R-parity conserv- 
ing event at the LHC in which supersymmetric particles were pair- 
produced. The colliding protons are shown coming in from the 
left and right. The collision has pair produced two massive susy 
particles, Ci o,nd (2 (dark blue). Each of these has been shown 
decaying to something visible (a or (3) and to an undetected neu- 
tralino (p or q). The typical event will also contain some initial- 
or final- state radiation, or other debris, represented here by g. In 
this figure it has been assumed that g consists entirely of visible 
particles. 



• the masses of their 'parent' particles are unknown, 

• the center-of-mass energy of the colhsion is not known, and 

• the boost along the beam axis of the collision center-of-mass is not 
known either. 

An example of such an event is shown schematically in Figure ^ where 
a pair of supersymmetric particle have been produced, each of which has 
decayed to some visible and some invisible daughters. 
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An important question to ask is "What model-independent information 
about sparticle masses can be deduced from events of this type?". The 
question can be seen to be a harder version of a number of older problems 
with which we are more familiar. 

Searches at the Large electron-positron collider (LEP) for pair produced 
sparticles have much in common with this problem, although to first order 
they did not suffer from the latter two of the above problems. In some ways, 
the problem at the LHC is more akin to that faced by the other hadronic 
collier experiments, such as UAl, UA2, CDF or DO where the Pi^-mass has 
been measured from its decay to a lepton and a neutrino. This has been 
achieved using the 'transverse mass' event variable, mr, a variable which on 
an event-by-event basis generates a lower bound for the W-mass, and the 
end-point of whose distribution is the W-mass ( [D El El ll| ) ■ Their case was 
easier than that at the LHC, however, as they only had one such decay per 
event, and in addition they could assume knowledge of the masses of both of 
the particles into which the W decayed. 

The approach to the problem first proposed in and subsequently devel- 
oped in jHj and [Zj proposes the creation of a new kinematic variable, mT2,^ 
analogous to the transverse mass, whose kinematic endpoint carries model 
independent information about (to first order) the mass difference between 
the primary and the secondary supersymmetric particles. Mention is made 
in [2] of generalisations to this variable (m^s, mr4, ...) which may be used 
when events contain extra missing particles (e.g. neutrinos) as well as the 
two neutralinos. 

The purpose of this article is not to discuss new physics results which 
might be obtained with mT2 (for these the reader is encouraged to read [Hj 
and [Zj) but rather it aims to take a closer look at more technical issues 
concerning the use and interpretation of mT2, and its related variables. It is 
hoped that by concentrating information on mT2 in this way, this article can 
act as a repository of mj'2 'know-how' for future investigations. 

2 A concrete example 

It is perhaps easiest to introduce and motivate the definition of the Cam- 
bridge mT2 variable using a concrete example. This allows the ingredients 
that make up ■mT2 to be introduced, one at a time in an almost 'natural' way. 
Readers who would prefer a 'top down' description of mT2, i-e. a description 



^Because of its use in supersymmetric events, mT2 has acquired the nickname of the 
'stransverse' mass. 
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which starts with a definition and then works towards its consequences, are 
directed to skip to section IHl where this approach is taken. 

The concrete example which will be used here is taken from ^ . This pa- 
per considered an (anomaly mediated) i?-parity conserving supersymmetric 
model whose key property was that it predicted a lightest chargino nearly 
mass degenerate with the lightest neutralino. With particular choices of 
model parameters, the only chargino decay mode available was: 



XI 
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Events containing two such decays, i.e. events containing two simultaneous 
decays of an unseen particle of unknown mass into another invisible particle 
of unknown mass and visible particle, are exactly the sort of events that we 
hope to analyse with 171x2- This we shall now begin to do. 

Considering for the moment just one of the decays of the form (^, one 
can write the Lorentz invariant statement 



ml+ =ml + w?o + 2 



El^E^^ cosh(Ar/) - ■ 



(2) 



where p^^ and p^^ indicate pion and neutralino 2-vectors in the transverse 
plane, and the transverse energies are defined by 



771 7r 
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p-TY + ml and = ^/(p^^)^ + m^o . (3) 



Also 



is the true rapidity, so that 



77 = - lo^ 



E + p, 



(4) 



tanht] = pz/E , sinhi] = Pz/Et , cosh t] = E / Et ■ 



(5) 



In a hadron collider, only the transverse components of a missing par- 
ticle's momentum can be inferred, so it is useful to define the transverse 
mass, 

m|(pj, p^"; m^o) = ml+ + m^o + 2{E^E^' - p}. ■ p^") (6) 

which, because cosh(x) > 1, is less than or equal to the mass of the lightest 
chargino, with equality only when the rapidity difference between the neu- 
tralino and the pion, Ar/^o^ is zero. All other Ar] lead to rriT < so z/we 
knew the neutralino momentum, we could use to give an event by event 
lower bound on the lightest chargino mass, was has been used this way 
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in the measurement of the mass. 

In events considered in this example, however, there are expected to be 
two unseen lightest supersymmetric particles (LSPs).^ Since only the sum of 
the missing transverse momentum of the two neutralinos is known, the best 
that can be done is to evaluate the quantity 



mm 



max{m|(p^''\^4^^^;m^o), m|(p^'", ^^J?^; m^o)}] (7) 



which is thus a lower bound on the square of the transverse mass, for 
events where two decays of the type occur. Note that this minimisation 
has forced us to introduce a pair of dummy two-vectors ^j:^ and which, 
constrained by the minimisation condition, parametrise our lack of knowledge 
about the two true neutralino momenta. Finally, we must recognise that 
under most circumstances, the value of m^o is unlikely to be known, or may 
only be known with limited precision. In order to make our ignorance of m^o 
explicit, we thus define a new free parameter, x, calling it the 'neutralino 
mass parameter', intending it to denote any guess we might have as to the 
true neutralino mass m^o. Using it in place of m^, we convert (jZj) into the 
following definition of a new kinematic variable: 

^T2{x) = mm [max|m|(p^''\^5^^^;x), ?«|(Pt'\ x)} 

T 



The quantity defined in (jH)) is the Cambridge 771^2 variable which is the 
subject of this document. 

Staying within the framework of this example, we can now go on to 
describe some of the the desirable model-independent properties which mT2 
possesses. 



2.1 Properties of mT2{x)- 

Firstly, is worth noting that the variable is not strictly a 'variable', and 
would more correctly be termed a 'function', as it retains a dependence on 
the unknown parameter x- Ideally, x would ideally be set equal to the mass 
of the missing heavy particle, but in most of the situations in which the 
variable is likely to be used, the mass of the invisible object is unlikely to be 
known, or may only be known with a large uncertainty. The x dependence 
remains, therefore. A more detailed discussion of how this can affects the 
use of mT2 takes place in section 12.2.21 

^Though there may also be other unseen particles - see section 12. 2. II 



5 



Secondly, from its method of construction, it is clear that for any given 
event 

rUj, + m^o < rnT2{fnxi^ — ^xt ^ 

m^ + X < "^T2(x) (10) 

It is certainly not immediately clear, however, that events can always exits 
for which mT2 is capable of reaching all of these endpoints. In fact it turns 
out that such events do always exist, and proof of this is given in section l!01 
So, having defined the quantity m™2^(x) by 

m??r(x) = max [mT2(x)] , (H) 

many events 

the important result to draw from all of this is that the upper kinematic limit 
of satisfies 

m-r(m^o) = m^+. (12) 
This is the main model-independent statement that mT2 is able to offer. 



2.2 Going beyond pairs of two body decays 

The scenario in which mT2 has been introduced, thus far, is relatively simple; 
each event contains a pair of charginos, and each of these decays via a two 
body decay into a charged pion and an unseen neutralino. We will now 
consider in more detail what happens when: 

• the neutralinos are not the only missing particles, 

• the initial (e.g. chargino) decays are not both two body decays, and 

• ^T2(x) is evaluated at values of x 7^ "^x?- 

2.2.1 Extra missing particles and multi-particle decays 

The need for mT2 to be adaptable to situations in which the neutralinos 
are not the only unobserved final-state particles may again be demonstrated 
using as an example the model of [7]. In this model, there were found to be 
some regions of parameter space in which three-body chargino decay, 

- l^^xl (13) 
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Invisible final state 



Figure 2: Schematic representation of a R-parity conserving event 
at the LHC in which supersymmetric particles are pair-produced. 
The colliding protons are shown conning in from the left and right. 
The collision has pair produced two massive susy particles, Ci 
and (2 (dark blue). Each of these has been shown decaying to a 
collection of visible particles (ai or (3i) and to a set of undetected 
particles (pi or qi). The purpose of the large spherical blobs is 
to hide the details of the decay process (es) involved; in principle 
they may contain anything, from one large n-body decay, to n — 1 
successive two-body decays. The typical event will also contain 
some initial- or final- state radiation, or other debris, represented 
here by g (the visible component) and g' (the invisible component). 
Comments in the text apply only principally to events in which g' 
is small enough to be neglected. 

had a rate comparable to that of the two-body decay (0) which we have 
aheady seen. Here, the presence of the neutrino (or antineutrino) in the final 
state means that we have even less information about the event than before. 
Nevertheless, one would like to benefit, if possible, from events in which one 
or two of these decays occur in place of the usual two-body decays. This 
type of event is just one of the general class of events depicted in figure El 
It is clear that one can immediately generalise the mT2 of (jH)) to suit 
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events like those in figure |21 in the following way. Define the new variable 
rriTx by: 



consistent splittings 



max < ai + ^ Pjf, (II A + II (Ij) 



■ (14) 



The phrase "consistent splittings", describing the constraint on the overall 
minimisation, needs a little explanation. There are two sets of unknown 
momenta. The first of these is JF = {p.i} U {qj}-, containing the unknown 
momenta of all the unobserved final-state particles. The other set, Ti. = {hi}, 
contains the momenta of any on-mass-shell particles which were present at 
an intermediate stage during the decays of the initial pair of sparticles (Ci 
or C2) to their final states. In other words, Ti. contains the momenta of any 
intermediate particles hidden within the large blobs in figure |21 Minimisation 
over "consistent splittings", then, means minimisation over all pi, qj G and 
all hi eH subject to: 

• all Pi, qi and hi being on their respective mass-shells, 

• momenta being conserved at all 'hidden' vertices in which a short lived 
intermediate particle with momentum /ij G 7i decays, and 

• the transverse components of B = J^iPi + J2j qj being consistent with 
the measured missing momentum pi^. 

It is because the last of these requirements that we need events in which g', 
the momentum carried by any invisible particles which are not descendants 
of a supersymmetric particle, (see figure ^ is negligible. Were there to be 
a large tail in the distribution of g', this would degrade the performance of 
rriTx and mT2- 



Example 

We illustrate the remarks of the previous section by returning to the example 
of [Zj in the case where where charginos could decay either by the two-body 
decay of ((H) or the three-body decay (fT!?jl . We can categorise events in 
this scenario by the number of missing particles in the event. When both 
charginos decay via we only have two missing particles (the neutralinos) . 
For each three-body decay which takes the place of one of these two-body 
decays we gain an extra missing particle in the form of a neutrino (or an- 
tineutrino). In short, the three categories of events could be summarised as 
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Figure 3: Simulations o/mTx("^x?) ~^xi f^''^ X = 2,3,4 using a 
simple phase-space Monte- Carlo generator program for a pair of 
decays q XiQ followed by Xi ~^ Xi ''^ or xt ~^ Xi ^ ^e- As the 
number of invisible particles increases, the proportion of events 
near the upper limit decreases. Within the figure, subscripts are 
indicated by square brackets. 



those containing one of the following: 

X^X^ - {vr^X?vr^X?, or e^uxWx'i, or e±z/x?e±^X?} ■ 

The events had been produced by a phase-space-only Monte-Carlo generator. 
Three distributions of the quantity rriTx, defined in dHj), were then generated 
from each of these sets of events. Using the number of missing particles to 
categorise these events, the values of ttitx measured in each case are referred 
to as mT2, "^Ts and mTi- The resulting distributions for m^xl^x?) ~ "^x? 
are shown in figure El 

It has already been mentioned that a key property of mT2 is that the 
kinematic endpoint of its distribution occurs at m™2^(m^o) = and so 
it is reassuring to see in figure El that a large number of events reach this 
endpoint. In the vicinity of the endpoint, the edge is seen to be sharp and 
near vertical. This shows that at the partonic level a measurement of m™|^ 
would provide an excellent constraint on the masses of the sparticles involved. 
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In section IHplots from [HllZl, which include realistic detector effects, will show 
that the subsequent smearing of the mT2 edge, while significant, is still small. 

Looking next at the mrs and ttita distributions, it is clear that the event 
fall-off in the vicinity of their kinematic endpoints is much less steep than 
in the case of mT2- This is hardly surprising, given the reduced amount of 
information available in these events. Later, in section 13.41 the relative frac- 
tion of events in the vicinity of the edge will be seen, more quantitatively, 
to be due to the larger number of simultaneous conditions that events near 
the edge must satisfy. Although the endpoint, itself, becomes increasingly 
harder to detect as the number of missing particles increases, the rriTx dis- 
tributions are all capable of inferring the mass scale associated with (in this 
case) m + = m^o from the overall widths of their distributions, which each 
scale with the endpoint position, albeit with some dependence on the decays 
themselves, and on factors such as the detector acceptance over the width of 
the distribution. 

Finally, one notes that the m^a distribution has sharp peak at = 
m^o + m^r, not seen in the mT2 and rriTA distributions. It will be shown 
later, in section 13. 5^ that this is an effect which can occur whenever the 
hypothesised decays on each side of the event are different.^ 

2.2.2 Other values of x 

Now we return to a brief look at the effect of evaluating rriTx distributions 
at values of x different to the true neutralino mass. 

Figure El shows the same data as in figure El but in addition it shows 
the distributions that would be obtained by evaluating rriTxix) values of 
X = fn^o ± 10%. In this particular example, where m^o = 162 GeV, 10% 
(16 GeV) errors in x result in similar fractional errors in AM^^ i.e. of a few 
tens of MeV. This shows that mT2 can be sensitive to small mass differences. 
In this example, too, we see a positive correlation between the change in 
X and the change in the position of the endpoint. These examples are not 
always typical, however. For example, in |3] the authors considered mT2 in 
the context of a double slepton decay to lepton a neutralino at SUGRA Point 
5, one of the five supergravity points proposed at [Hj and described in pj. In 
this model, the difference in mass between the decaying and final sparticles 



■^'Different', in this context, means 'being such that the minimum total invariant mass 
attainable by the particles on one side of the event is not equal to the minimum total 
invariant mass attainable by the particles on the other side of the event'. This happens 
principally when the particle content of each decay differs. In the case of tot3 in, the 
AMSB example scenario, the two dissimilar minima are + m^o and me + rrii, + m^o . 
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Figure 4: The distortion of rrixxix) " X when the LSP mass 
parameter, x> varied by ± 10% about the 'ideal' value of 
m^o. These curves show that rnxxix) ~ X remains sensitive to 
the mass difference AM^, = m + — m^o. In this simulation 



AMj^^ = 0.845 GeV, m^o = 161.6 GeV, and the electron and 
neutrino mass were neglected. The normalisation is arbitrary. 
Within the figure, subscripts are indicated by square brackets. 



(157.1 — 121.5 = 35.6 GeV) is approximately 40 times larger than in the 
AMSB case, and so at SUGRA Point 5 we see a negative correlation between 
changes in x and mT2{x) ~ X- This is illustrated in figure El Differing kinds 
of behaviour, such as these, are typical of a variable like 771^2 which has input 
scales (e.g. m^r, ^^+, X " ) which can have a large number of 

relative hierarchies associated with them. For example, the AMSB scenario 
has ~ m+ — X ^ X ^ '"^v+ ; while SUGRA Point 5 has <^ m,- — x ^ 

Al Al 

X < mj. 

We now take a final look at how depends upon x by looking not just 
at events near the kinematic endpoint, but at events in general. To help, we 
define y{x)^ ^ rescaling of mT2(x)) as follows: 



By looking back at Q and (fT^ . this variable can be seen to map mT2{'m-x^) 



y{x) 



mT2{x) -X-m^ 



(15) 




'TT 
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25 50 75 100 125 150 175 200 

m, (GeV) 

Figure 5: Variation of with x for a set of l^l^ —>■ 

l~^Xi^~Xi events generated by a phase-space Monte- Carlo using 
mi = 157.1 GeV and m^o = 121.5 GeV. Note that m?;|^(x) - x 
decreases as x increases. The asymptote has unit gradient. 



into the range [0, l].'^ This makes it easier to compare values of mT2 coming 
from events with different sparticle masses. A value of ?/(m^o) close to 
(or 1) indicates an event close to the lower (or upper) kinematic endpoint 
of the mT2ijn^^i) distribution. Figure IHl shows how yix) varies with x for 
ten random events, each generated using a random set of masses satisfying 
m,r + fTi^o < as described in the figure caption. The main conclusion 
to draw from these plots is that there is no easy way to say, in advance, 
how rnT2{x) will vary with x given event, even in the vicinity of m^o. 
In general mT2{x) can rise, fall, or even be stationary w.r.t to x near m^o, 
depending on the masses of the particles involved in the decays. 



^Note that 771^2 (x) is in general not constrained to lie in this range. 
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Figure 6: These plots show examples of how y{x)j defined in / f73|) . 
can depend on x- The plots were generated using the following 
procedure. Ten sets of masses satisfying m.^ + m^o < were 
randomly generated. According to each set of masses, a phase- 
space Monte- Carlo generated a single event of the type shown in 
figureUl containing two xf ^^X? decays. The plots above show 
how, in each event, the value of y{x) (a dimensionless rescaling 
of mT2{x)) depended upon x over the range < x < 2m^o. The 
true value of the neutralino mass (i.e. that used in the Monte- 
Carlo for the decay) is marked by the vertical line at the center of 
each plot (cyan), while the other vertical line marks the value of 
the pion mass (magenta). The short and long vertical ticks (dark 
blue) mark m + /2 and m + respectively. 



3 Some mathematical results concerning the 
variable rriTx 

In section mT2 and its friends were defined using lab-frame momenta and 
with minimisation conditions (such as + = not specified in 
lorentz-invariant forms. The definition of from section |21 may be sum- 
marised as follows: 

"^T2(Pr,Pr)f^r;x)= min [max {m^(p2i, fi^; x), m^(p^, 1^2; x)}] , (16) 



13 



where 



mUPT, Pl; X) = ml + + 2{ElrE^ - ■ p|), (17) 



in which Elr 



Pt 



2 



I 5 



Pt 



+ X^- III the following sections we 



will replace these definitions by equivalent, but explicitly lorentz- invariant 
ones, which are easier to manipulate mathematically. 

3.1 Definitions 

The natural way to write in a manifestly (1, 2)-lorentz invariant form is 
as follows: 

771^2(0^5 ^5 x) = iniii max {(a + pY , {P + qY} 

( p + q = yiA \ 

I P^ = q = X^ J 

Here, 771x2 has been written as a function of the four 1 + 2-dimensional 
lorentz vectors which describe each event (a, /?, S and A) and one real pa- 
rameter X- The transverse lorentz vectors of the two visible particles coming 
from each of the hidden decays are represented by a and (3, while S represents 
the total transverse energy- momentum seen in the event. This is consistent 
with the notation used in Figure ^ The only new vector. A, defines the 
laboratory frame by being the (1, 2)-energy-momentum of a particle of unit 
mass (A^ = 1) at rest in the laboratory. The total transverse momentum 
of the event (visible and invisible) can only be assumed to be zero in the 
laboratory frame, and so knowledge of how to boost to the laboratory frame 
is essential. This is why A is needed. 

Note that p8|) includes a minimisation over ^/s, a parameter which ac- 
counts for our lack of knowledge of the center-of-mass energy of the whole 
event. The requirement that the hypothesised neutralino momenta are real, 
i.e. (y^A — S)^ > (2%)^, constrains y/s to be chosen from the region in which 



v^> A.S + y'4x2 + ((A.S)2-S2). (19) 
Similarly, one can also define m-rx in a manifestly lorenz-invariant form: 



m^j^(a, (3, S, A; x) = niin max {(d + pY, {(3 + g)^} 



over 



p + q = v/iA - S, 
Pi^qi and hi all on their mass shells, and ^, (20) 
momenta conserved at all internal h decays 
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in which the same notation has been used as in figure El and equation (fT^ . 
and in which the 'hats' indicate summation over all vectors of a set (e.g. 

If desired, one may remove the 'max' at the expense of moving to lorentz 
four-vectors and remembering to minimise over all possible longitudinal boost 
of the center of momentum (here denoted by the lotentx boost L^): 



This way of representing mT2 niost clearly captures the spirit in which it 
provides an event-by-event lower bound on the initial sparticle mass. 

3.2 Results concerning mT2 

In the case where both visible particles have the same mass, i.e. in the case 
where = = mf, (|T8|l may be re- written in the form: 



It was shown in [7j that the solution of must select vectors p and 
q for which a.p = j3.q. Using this information, one may perform half of 
the minimisation in ()24|) analytically. This allows ()24|) to be re-written as a 
minimisation over a single real variable, ^/s, as follows: 




where (21) 



(22) 
(23) 




2max{a.j9,/?.g} .(24) 



m%^{a, (3, S, A;x) = m^ + + min {(x.B)Q - 



s 



a^Q - Amj^B^Q - 4^^ , 



where 



Q 



1 - 



(A.5)2 




(26) 




a + /3, 
a — (3, 
v^A-S. 



(27) 



A 



and 




B 



(29) 



15 



We should note that the constraint which has just been imposed, namely 
a.p = f3.q, is more stringent than the ^/s constraint (|TIHl which was only 
there to ensure that the hypothesised neutralinos were not tachyonic. As a 
consequence, the range over which ^/s may be varied when performing the 
minimisation in ()25p must be replaced by the stronger condition that each 
of the quantities under radicals in be positive. 

It is interesting to note that if we define two new transverse lorentz vectors 
{a and B) via a rescaling of existing transverse lorentz vectors according to 



a = cryQ, and (30) 
B = B^, (31) 

then we can rewrite in the form 

1 



771^2 (a, P, S, A; x) = m] + + 7; mm 



a.B) - Ja^ - Amf^B^ - 4% 



.(32) 



This is not much of an improvement in itself, but it motivates the definition 
of two new lorentz four-vectors; 



a = {(J , J a'^ — Amf) , and (33) 



B = {B,^b'-Ax'), (34) 

which we see, by construction, satisfy the following fixed-mass relations: 

= \^ = 2m;, and (35) 
= = 2x. (36) 

In terms of these new lorentz four-vectors, then, we can finally re-write (jH^ 
as 

m|2'(a>/^>^.^;x) = "^f + + ^ min (o-.e) , or (37) 

mT2[a,i5,"L,k]x) = ]-m.m a + B . (38) 

2 

It is interesting to note that the constant mass relations ()35j) and ()36j) . taken 
together with the definition of 771x2' shown in (jHH|) . make it self evident that 
the value of 771T2' obtained in a given event is bounded below by m/ + x, as 
expected. 
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Approximations 



To get a better idea of the way in wiiicli 771^2 depends on its inputs, one 
might hope to find a concise closed-form analytic definition of the variable. 
Thus far, however, mT2 and mT2 have resisted all attempts to write them 
in forms simpler than (|T6|) . (|T8|l and (j25|l . except in a few special cases. For 
example, in the special case of events in which the spatial part of the total 
visible transverse momentum is seen to be zero in the laboratory frame (i.e. 
events for which (S.A)^ = S^) one can show that ()25j) is equivalent to: 



ml," {a, (3, A; x) ^ + + xv'4(a.A)(/?.A) - (-A^) (39) 



in the laboratory frame ). 

The limit of validity of ()39|1 can be explored as follows. The laboratory 
frame energy that this special case solution assigns to p and q is given by 

P.A = «.A = , (40) 
^4(a,A)(/J,A) - 

and so the velocity of the boost needed to take the laboratory frame to the 
one in which the invisible particles are back to back could be written, in this 
special 

(0 = ) |v|^ = pl/ipA + qAf (41) 

{(S.A)^-Sn{4(a.A)(/3.A)-(-A2)} 



Aia.Ayx' 

pUEgEp + m^ + p^.p^T) 
2{E^ + EpYx^ 

in the laboratory frame ) 



(42) 
(43) 



In the light of the above, we can interpret ()39|) as the leading term in an 
expansion of in powers of |v|^ as defined in ()42|). Given a particular 
event, all the quantities in (jl^ may be evaluated, so one can safely use (jHIIjl 
to evaluate m^, events in which |v|^ is observed to satisfy |v|^ <C 1.^ 



^The reader is warned not to mistake |v| for the speed associated with an actual boost 
(real or conjectured) connected with the neutralino pair; |v| could, for example, even 
exceed the speed of light if x were made sufficiently small! It should only be assumed 
that as |v| — > 0, |v| will tend to the speed, in the laboratory frame, associated with the 
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X 



X 



X 



(a) (b) (c) 

Figure 7: A diagram demonstrating that the minimisation over 
some parameter of the maximum of two well-behaved functions 
may occur either at (a) a minimum value of one of them, or (b) 
when they are equal, or (c) at the boundary of the domain. 



3.3 Extremal values of mT2 

In this section, we show that the maximum value which mT2{'n^x'i) ^^"^ attain, 
for a given set of particle masses, is indeed the mass of the initial sparticle.® 
We start from definition (jHl). We also describe the region of decay-phase- 
space which contains events which occur close to this kinematic endpoint. 

To find the range of values mT2 inay take we first let /i = m|,(p J*^* , ; m^o ) , 

and /2 = m|^(py^^\ m^o). We then note that the minimum over a pa- 
rameter X of the maximum of fi{x) and f2{x) can occur at a local minimum, 
f[(2){.x*) = 0, provided fi(2){x*) > f2{i){x*), as shown in figure IH7^ . Alter- 
natively the minimum can occur when the functions cross one another when 
fl = f2 (figure ElSb) or at a boundary (figure ElSfc) • The parameter x corre- 
sponds to the fraction of the the missing momentum (in one of the transverse 
directions) which is assigned to each half of the event. Since /i, /2 — ^ c)o as 
X —>■ ±oo figure 13.3b is not relevant to our minimisation problem. 

To find which of (a) or (b) is pertinent, consider an unconstrained min- 



energy-momentum vector B. (This is the B which was originally defined in (|29|l and whose 
y/s value was selected by the minimisation process in (1^ .^ 

^Up to this point, within the context of the AMSB example, it has only been shown 
that 171x2 ) bounded above by ^1^+ . It has not yet been shown that 771^2 can attain 
this bound. The purpose of this section is to show that it can. 
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imisation over of mf{p^, m^o). Using the relationship 
where Ip"^ = ^ + ?77. 0, it is straightforward to show that, 



(44) 



9m: 



2\E1^^-Pl\ k = l,2. (45) 

This means that at an unconstrained minimum of rn^ we have 

V? = yir , (46) 

where we introduce the notation v^- = pt/Et, = ^t/^t, in which p^- 
and represent the true transverse momentum and velocity of a particle, 
while and yir are assigned by the minimisation. 

Using the basis (t, x, y) with the metric diag(l,-l,-l), one can write 

m^ = (E^°\pr)-(E^?\p^°*) , (47) 

where E^^ = E^ + and p^* = Pt + ^t- This 1+2 dimensional Lorentz 
invariant can be evaluated in any frame boosted from the lab in the transverse 
plane. PU)) has told us that at the unconstrained minimum the transverse 
velocities and yt^ are equal; a statement necessarily true in all transverse 
frames, including the special one in which both the transverse velocities (and 
associated momenta) are zero. Evaluating ()47|1 in this frame, we find that 
the unconstrained minimum of (jTfjl then becomes (m,r + ^x?' ' (^^^ 

m^o, 0, 0), and we recover the expected result 

m^i" = m^ + m^o . (48) 

We therefore conclude that the function has only one stationary value 
and it is the global minimum, and is common to both sides of the event 
provided the same type of particles are emitted. Thus when /i is minimum 
it cannot be greater than /2, and so the minimisation in (jHJ forces /i = /2- 
This could of course occur when both fi and /2 are at their global minima, 
in which case mT2 takes its minimum value: 

mT2™''' = m^ + m^o . (49) 

To summarise, when the same particles are emitted from both sides of 
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the event, mT2 may be defined as the minimum of m>p subject to the two 
constraints m^-* = and i^^"* + i^^"* = ^j,. The condition for the minimi- 
sation can be calculated by lagrange-multiplier methods, the result of which 
is that the velocity vectors ^i^'^"* of the assigned neutralino momenta ^^''^^ 
must satisfy 

(>i«-vf') cx i^P-^rf). (50) 

To find the maximum of mT2 over many events we note that for each 
event the minimisation will select hypothesised momenta satisfying (jSUj) . We 
now note events can occur in which the true transverse velocities of the 
neutralinos were exactly those which were assigned by the minimisation, i.e. 
they can satisfy 

vjS'" = 4", v*' = 0<?'. (51) 

These events will have both hypothesised transverse masses equal not only 
to each other but also to true transverse masses which would have been 
calculated if the neutralino momenta had been known: 

m?(p^",fi^°«)=m«(pf ,#) (52) 

If events occur where, in addition to the transverse components of the neu- 
tralino momenta satisfying (jM^ . the rapidity differences satisfy ^;^o(i) = 77^(1) 
and ^7^0 (2) = ?77r(2)5 then by ^ mT2 will equal the true mass of the chargino. 
Combining this with ()49j) and recalling that mT2 cannot be greater than the 
chargino mass by construction, we can see that the event-by event distribu- 
tion of mT2 can span the range 

m^o-Fm^ < mT2 < ^xt' ^^'^^ 
3.4 Extremal values of ttitx 

In the last section we looked at the conditions under which events can gen- 
erate mT2 values near the kinematic endpoint. Here we will look at some of 
the ways these conditions become modified for m^s and rriTi events. 

Consider once again events from the AMSB scenario in which a chargino 
is produced and then decays to x\ ^ ^e- If we expand the Lorentz invariant 

{m^+f = {p^o+Pe+p,f (54) 

we obtain three mass-squared terms for each of the decay particles and three 
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cross-terms. The cross-terms can each be written in the form 



2pa ■pt = 2 [e^t^^E^^^ cosh(Ar7,b) - ■ p?] , (55) 

hke the cross term in Q. If the neutrahno and neutrino transverse momenta 
were individually known we could evaluate the transverse mass, 



{E^E^ - Pt ■ Pt) + {EtE^ - Pt ■ Pt) + (^t^t - Pt " Pt) J , 

(56) 

where the neutrino mass is assumed to be negligible. rriT will be equal to the 
xt niass in events where Arjab = for all pairs of e, z/g, and Xi- 

Using, in (jHl), the three-particle definition of rriT from instead of the 
two-particle definition ©, one defines mT4, the analogue of for the case 
of four missing particles. The constraint on the unobserved momenta will, of 
course, have to be modified to read 

q^^) + q^(^^ + + q^^^^ = , (57) 

where the labels (1) and (2) indicate which chargino the particles were emit- 
ted from. 

The conditions for the minimisation required to calculate can be 
calculated just as for mT2- The Euler-Lagrange (E-L) equations involving 

^ a.d ^ (58) 

show that the minimisation will select the invisible particles' momenta such 
that The other E-L equations reproduce (|Hn|) but with elec- 

trons replacing pions. 

This means that when calculating one can replace the missing par- 
ticles from each chargino decay with a pseudo-particle with mass equal to 
the sum of the masses of those invisible particles and proceed as for mT2- In 
the case of leptonic chargino decay the mass of the neutrino can be safely 
neglected in comparison to that of the Xi? and the constraint u^*-*'' = u^*-*^ is 
equivalent to q^*-*^ = (0,0). 

The distribution over events of rriTA will have fewer entries near the upper 
kinematic limit (mT4 = fT^^+) because when more particles go undetected an 
event at that limit must satisfy a larger number of constraints. For fully 
leptonic chargino-pair decay, there are six constraints of the type At] = 0, 
two p'^^^ = and finally the modified constraint from (jSUjl . This effect can 
be seen in figure El for events where a total of two, three and four invisible 
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particles are produced. 



3.5 Asymmetric decays 

In the preceding two sections we have seen that when the decays on each side 
of the event are the same (i.e. both initial sparticles decay to the same set 
of daughter particles) then the resulting kinematic variables, 171x2 and rriTA, 
have very similar properties. The only significant difference we have seen 
is the reduced density of events near the upper kinematic endpoint of 771^4 
relative to mT2- Why, then, is the rriTs distribution, shown in figure El seen 
to have a shape significantly different from the and 771^4 distributions? 
Specifically, why does it have the strong peak at low values not shared by 
the other two? 

The difference occurs because the visible particles on each side of an m^s 
event are different (on one side y^^e^v and on the other to x^^t^^) and so 
the unconstrained minima of the values of tut on each side of the event are 
not equal as they are in the case of mT2 and uita- 



It is thus possible for some of the events can then fall into the category shown 
in figure ISiSK, producing a peak of events with ttitx = f^x" ~^ 

4 LHC case studies using 771^2 

In this section we highlight some physics studies for the LHC which demon- 
strate that the background and the detector effects do not prevent mT2 from 
being a useful experimental variable. We investigate points from three dif- 
ferent models, under two different classes of mass hierarchy. 



The first two points discussed are the mSUGRA point 5 (S5) and a point from 
the optimised string model (01) discussed in [Hj. The relevant parameters of 
these models are 




(59) 



4.1 Case 1 - mSUGRA-like points 



^3/2 

m^/2 



300 GeV, mo = 100 GeV, Aq = 300 GeV, tan (3 = 2.1, 
250 GeV, tan/? = 10, ^ = 7r/4 
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AM (GeV) 



(a) S5, hard cuts 




AM (GeV) 

(b) 01, hard cuts 



Figure 8: AM distributions obtained at S5 and 01 after applying 
the cuts to a 100 fb~^ sample of crCr, fiRfin, e^ei and jlLp'i 
dislepton events. Events from light right-sleptons (unhatched) are 
stacked on top of those from heavier left-sleptons (hatched). Only 
events from OSSF leptons combinations are shown. The plots are 
generated without OSDF background subtraction, but were it to 
be performed, no significant differences would be apparent as only 
4 (12) signal events are able to pass OSDF soft cuts at S5 (01). 
Arrows indicate the values of AM™'^^ predicted by theory for the 
two types of slepton in each model. A red vertical line is drawn 
through each plot at half the W mass. 



respectively and /x > in both cases. For these points we are looking at 
dislepton production from a hard process which decays as — > xV'^ i ^-iid so 
the mass-hierarchy is <^ m^- — x ^ x < 

For S5 and 01, all events, except the qq W^W~ background processes, 
were simulated by HERWIG-6.0 JT^. The W-pair events were generated by 
ISAJET-7.42 The events for these two points were generated at 100 

fb~^. This is expected to correspond to running at high luminosity for one 
year. 

Since there are two different processes being analysed, there are different 
cuts to apply. As this is not intended to introduce new physics, here we 
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present only the major cuts used. For more detail about the cuts and the 
techniques used, see jEllZl- 

The events used for S5 and 01 are required to have one opposite sign 
same family (OSSF) pair of isolated leptons with p^j. > 50 GeV and > 
30 GeV. These events cannot contain any other isolated leptons. Also, events 
containing one ore more jets with pip > 40 GeV are vetoed. This helps reduce 
the standard model backgrounds. 

The variable AM is defined as 

iAM)^^{MUmi)y~ml (60) 

This variable is what is studied for the points S5 and 01, for reasons given in 
in]. The desired dislepton events have very little jet activity and the dislepton 
production cross sections are typically two orders of magnitude smaller than 
the squark/gluino production cross sections. There are also irreducible SM 
backgrounds (primarily W~^W~ — > l^l'vv and ti —>■ bbW^W~ —>■ 
in cases where jets are below the reconstruction threshold or are outside 
detector acceptance) which have signatures identical to dislepton events. The 
smallness of the signal and the presence of these backgrounds would cause 
problems for naive straight-line fitting technique. Instead, the technique 
described in [B] is used for the estimation of the edge precision. 

Figure |H1 shows the AM distributions obtained at S5 and 01 after ap- 
plying the cuts to a 100 fb" sample of signal dislepton events. Events from 
lighter sleptons {cr and /ii?) occupy the unhatched region in each plot, while 
the events from heavier left-sleptons are cross hatched. It will be noted that 
events from both light and heavy sleptons succeed in passing the cuts in both 
models. In principle, then, there are two edges to be observed in each of the 
models: one for the lighter slepton and one for the heavier. We note, however, 
that as the slepton masses increase, their production is strongly suppressed, 
and so there are very few heavy slepton events at 01 where there are in fact 
none within 10 GeV of the kinematic limit. It is readily observed that at 
the three remaining edges, where statistics are higher, there is good agree- 
ment between the theoretical prediction and the observed endpoint of each 
distribution. 

Significant numbers of SM background events also pass the cuts. These 
can be well modelled by its opposite sign different family (OSDF) coun- 
terpart. As the signals from dislepton pair production are expected to be 
purely OSSF, we can use OSDF background subtraction. Supersymmetric 
backgrounds also have to be considered. Again, in this case the OSSF distri- 
butions are well modelled by OSDF events passing the same cuts. So again. 
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Figure 9: Flavour- subtracted AM distributions for combined sig- 
nal and background at S5 and 01 after applying the cuts. A red 
vertical line is drawn through each plot at half the W mass. Plots 
are shaded to the left of this line in order to draw attention to the 
events which reconstruct above this point. Compare these plots 
with those of Figure{^ which contain only signal events. The ar- 
rows of Figure are provided for comparison. 



different-family background subtraction is used. 

All events for S5 and 01 (signals and backgrounds) are combined in Fig- 
ureiniafter different-family background subtraction. The reader is encouraged 
to compare these plots with those from Figure |H1 showing the desired signal 
shapes. As expected, all signal shape information is lost to the left of mw/'2 
due to obhteration by the SM backgrounds. To the right of this point, at 
least one clear edge is observable in both models (the left-slepton edge at S5, 
and the right-slepton edge at 01) and in both sets of cuts. The hard cuts are 
able to suppress the supersymmetric backgrounds to such a degree that there 
is even compelling evidence at 01 for the existence of two edges, although 
the lack of statistics in the higher edge limits the precision with which the 
endpoint may be located. 
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Point 


(GeV) 


AM,, 
(MeV) 




Xl ^ XV^'^M 


SPS-300 


165 


886 


17.0 % 


15.9 % 


A-250 


101 


766 


15.4 % 


13.9 % 


SPS-250 


159 


1798 


21.9 % 


21.5 % 


A-200 


97 


1603 


22.5 % 


22.2 % 



Table 1: The lightest chargino mass, the mass difference, AM,, = 
—m^o, and two chargino branching ratios for the AMSB-like 
points discussed in section \4.S\ The hadronic branching ratios 
can be found in J^. 



4.2 Case 2 - AMSB-like scenarios 

The characteristic signature for anomaly-mediated supersymmetry break- 
ing is the near mass- degeneracy of the lightest chargino and the lightest 
neutralino. The xt therefore decays to a neutralino plus (relatively) light 
standard-model particles. For a small mass difference, AM,, = — m^o, 
the largest xt branching ratios are to to Xi^^^h where I E e,fi. 

The mass hierarchy, 

m^r or {mi + m^) ^ — x ^ X ^ ' 

is therefore very different to the previous case study. 

HERWIG-6.3 was used to generate 30 fb~^ of unweighted inclusive super- 
symmetry events. HERWIG was also used to generate the background. For all 
the points, the results were passed through the ATLAS fast detector simu- 
lator, ATLFAST^^l- The signal-enhancing cuts require missing transverse en- 
ergy, = 500 GeV, leading jet transverse momentum, p™},) = 400 GeV 
and transverse sphericity, S"™™ = 0.05. There are also cuts on the tracks, 
these are described in more detail in |7j. 

We consider AMSB-like points, which have the following parameters: 

mo = 450 GeV, mg/s = 60 TeV, tan /3 = 10, > 0, 
mo = 500 GeV, mg/s = 36 TeV, tan /? = 10, /i > 0. 

and for which the /i parameter has been adjusted at the electroweak scale 
in order to investigate different values of AM,,, as discussed in j7]. Some 
masses and branching ratios can be found in tabled 
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(a) SPS-300 (b) A-250 

Figure 10: The mT2 — f^xl distribution for (a) the point SPS- 
300, and (b) the point A-250. The signal consists of the two 
solid regions labelled + in the legend. The upper kinematic 
limit of mT2 ~^x° f^''^ signal events is marked with a dotted line. 
Note the sharp fall- off in the distribution near the kinematic edge 
at mT2 — = AMyj . 



The first two points have large branching ratios for the decay xf — > 'rc'^Xi- 
This means that chargino-pair decay can easily generate the topology shown 
in section ^ We therefore plot distributions of mT2 — ' -^"-"^ which signal 
events in a perfect detector would lie in the range [m,r, AM^J. The results 
(see figure irUj) show that could be used to measure the small mass 
difference between the Xi ^ind the Xi this model, provided the signal 
cross-section is sufficiently large. 

The second pair of points each have a larger leptonic branching ratio, and 
so for these points the fully leptonic channel was investigated. Since there 
are now four missing particles in the final state, (two neutralinos, and two 
neutrinos), distributions of rriTi — m^o were plotted. For a perfect detector, 
these are restricted to lie in the range [rrie/^, AM^-^]. 

The signal events are again indicated by the solid shades in the histograms 
in figure ^3 Again, it can be observed that the distribution lies within the 
expected range. The distribution is skewed to lower values because more 
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Figure 11: The mT2 — f^xl distribution for (a) the point SPS- 
250, and (b) the point A-200. The signal consists of the two 
solid regions labelled + in the legend. The upper kinematic 
limit of mT2 ~'^x\ ^^'^ signal events is marked with a dotted line. 
Note the fall-off in the distribution near the kinematic edge at 



particles are missing, and so more constraints must be satisfied for an event 
to approacli tlie upper limit (as was seen in figure E)). 

The sensitivity of rriTx to the estimated mass of the neutralino was shown 
in figure lU It has been found that rriTx shows similar insensitivity to mea- 
surement uncertainties in the missing transverse momentum vector. This 
behaviour can be (at least partially) understood from the non-relativistic 
limit of mT2, when the proportionality in ()50|) becomes an equality and 

1 2 

m|2-("^7r+m^o)^ = (m^i^T - "^x?Pt - "^x?Pt ) +C ((vt ■ Vt)^) 

4r /TZ77- 7TZ , 

(61) 

The low sensitivity to the (possibly poorly-measured) quantities m^o and fir 
comes from the fact that in (jF)T|l they are multiplied by the quantities pj" and 
respectively, which are both small in this mass regime. 
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5 Conclusion 



This paper has attempted to achieve three objectives. Firstly it seeks to 
introduce a new set of kinematic variables {mT2, ''^ts, •••}; which are specially 
designed to extract information from a particular class of troublesome events 
that we are likely to see at next generation hadron and lepton colliders. These 
events are those containing a pair of particles of identical (but unknown) mass 
which subsequently decay into groups of particles, each containing one or 
more invisible (possibly massive) particles. An example of this kind of event 
might be pair production of sleptons at the LHC, followed by subsequent 
sleptonic decay to leptons and neutralinos. Secondly this paper attempts 
to get to the bottom of these new variables; it describes the regimes with 
in which they can or cannot be trusted, develops useful approximations to 
them, and shows generally how one could go about calculating this variable 
for real. The approximations to the variables are not only useful in their own 
right, but are even more useful as guides which illustrate the dependence of 
the variables upon its inputs. Finally, this article seeks to show with a couple 
of examples, real use of these variables in physics analyses. These hopefully 
show that mT2 and its chums are able to provide vital and new information 
about particle masses from events that would at first glance appear to contain 
so many unknown quantities as to be useless. 

In conclusion, we believe that 771^2 is invaluable tool for physicists working 
at the LHC, and other future colliders, and we hope that this document will 
encourage its use. 
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